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ANALYSIS OF COMPOSITE PANEL-STIFFENER DEBONDING 
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1 .2 

Ronald Krueger and Pierre J. Minguet 


ABSTRACT 

Interlaminar fracture mechanics has proven useful for characterizing the onset of 
delaminations in composites and has been used with limited success primarily to 
investigate onset in fracture toughness specimens and laboratory size coupon type 
specimens. Future acceptance of the methodology by industry and certification 
authorities however, requires the successful demonstration of the methodology on 
structural level. For this purpose a panel was selected that is reinforced with stringers. 
Shear loading causes the panel to buckle and the resulting out-of-plane deformations 
initiate skin/stringer separation at the location of an embedded defect. For finite element 
analysis, the panel and surrounding load fixture were modeled with shell elements. A 
small section of the stringer foot and the panel in the vicinity of the embedded defect 
were modeled with a local 3D solid model. A failure index was calculated by correlating 
computed mixed-mode strain energy release rates with the mixed-mode failure criterion 
of the graphite/epoxy material. Computed failure indices were in good agreement with 
results from models where the entire delaminated section of the stiffener foot had been 
modeled with solid elements. The shell-to-solid connection influenced the computed 
failure indices and local refinement of the shell model across the stringer foot and web 
was required to improve the results. The study confirmed that the section modeled 
locally did not have to include the entire delaminated section. The use of a smaller local 
inserts reduced model size without compromising the computed failure indices. 


1. BACKGROUND 

Many composite components in aerospace structures are made of flat or curved panels 
with co-cured or adhesively bonded frames and stiffeners. Recent studies focused on the 
investigation of the debonding mechanism and included testing of composite skin/stiffener 
panels and failure analysis using shell models [1-4]. Over the last decade a consistent step-wise 
approach has been developed which uses experiments to detect the failure mechanism, 
computational stress analysis to determine the location of first matrix cracking and 
computational fracture mechanics to investigate the potential for delamination growth. Testing of 
skin gage stiffened panels designed for pressurized aircraft fuselage has shown that bond failure 
at the tip of the frame flange is an important and very likely failure mode [5], Comparatively 
simple specimens consisting of a stringer flange bonded onto a skin have been developed to 
study skin/stiffener debonding [6-8], The failure mode that initiates at the tip of the flange in 
these specimens is nearly identical to the failure mode observed in the full-scale panels and the 
frame pull-off specimens [7, 9, 10]. A methodology based on fracture mechanics [11] has proven 
useful for characterizing the onset and growth of delaminations in composites and has been used 
with limited success to investigate delamination onset and debonding in spimple laboratory 
coupon type specimens [9, 10]. Future acceptance of a fracture mechanics methodology by 
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industry and certification authorities however, requires the successful demonstration of the 
methodology on structural level. 

For the demonstration of a fracture mechanics methodology on structural level, a stringer 
stiffened panel as shown in Figure 1, had been selected and analyzed in a related study [12]. The 
square panel is made of IM7/8552 carbon/epoxy tape and reinforced with three stringers made of 
IM7/8552 carbon/epoxy plain weave fabric. The stiffened panel is bolted to a steel picture frame 
and subjected to pure in-plane shear loading. A finite element model is shown in Figure la. 
During manufacturing an artificial defect had been placed at the termination of the center 
stiffener. Sufficient shear loading causes the panel to buckle as shown in Figure la. The resulting 
out-of-plane deformation causes skin/stringer separation at the location of the initial defect as 
shown in the detail of the deformed mesh in Figure lb. A total of eight delamination lengths 
between 81.9 and 355.6 mm were modeled using short and long insert models as shown in 
Figures lb and c. The mixed-mode strain energy release rates were calculated using the virtual 
crack closure technique across the width of the stringer foot. A failure index was calculated by 
correlating the results with the mixed-mode failure criterion of the graphite/epoxy material [12]. 

The objective of this research is to establish that the section modeled locally did not have 
to include the entire delaminated section as shown in Figures 1 b and c. The goal was to use 
smaller local inserts and thereby reduce model size and CPU time without compromising the 
computed failure indices. A total of seven delamination lengths between 81.9 and 355.6 mm 
were modeled and computed failure indices were compared to previous results where the entire 
delaminated section had been modeled with solid elements [12]. 


2. METHODOLOGY 

2,1. Interlaminar Fracture Mechanics 

Interlaminar fracture mechanics has proven useful for characterizing the onset and 
growth of delaminations [11, 13]. When using fracture mechanics the total strain energy release 
rate, Gt, the mode I component due to interlaminar tension, Gi, the mode II component due to 
interlaminar sliding shear, Gn, and the mode III component, Gm, due to interlaminar scissoring 
shear, as shown in Figure 2, are calculated along the delamination. The calculated G/, Gn, and 
Gm components are then compared to interlaminar fracture toughness values in order to predict 
delamination onset or growth. The interlaminar fracture toughness values are determined 
experimentally over a range of mode mixities from pure mode I loading to pure mode II loading 
[14-16]. 

A quasi static mixed-mode fracture criterion is determined by plotting the interlaminar 
fracture toughness, G c , versus the mixed-mode ratio, Gn/Gj. The fracture toughness data is 
generated experimentally using pure Mode I (G///G7K)) Double Cantilever Beam (DCB), pure 
Mode II (Gh/Gt=\ ) four point End Notched Flexure (4ENF), and Mixed Mode Bending (MMB) 
tests of varying ratios as shown in Figure 3 for IM7/8852 carbon epoxy material. A failure 
criterion was suggested by Benzeggah and Kenane [17] using a simple mathematical relationship 
between G c and Gu/Gt 


G c = G Ic + { G IIc~ G Ic)' 



( 1 ) 
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In this expression G/ c and Gu c are the experimentally-determined fracture toughness data for 
mode I and II as shown in Figure 3. The factor rj was determined by a curve fit using the 
Levenberg-Marquardt algorithm in KaleidaGraph™. Failure is expected when, for a given mixed 
mode ratio Gu/Gt, the calculated total energy release rate, GY, exceeds the interlaminar fracture 
toughness, G c . Although several specimens have been suggested for the measurement of the 
mode III interlaminar fracture toughness property [18-20], an interaction criterion incorporating 
the scissoring shear, however, has not yet been established. The edge-cracked torsion test (ECT) 
is being considered for standardization [21, 22]. 


2.2. Analysis Tools 

2.2.1 Virtual Crack Closure Technique 

A variety of methods are used in the literature to compute the strain energy release rate 
based on results obtained from finite element analysis. For delaminations in laminated composite 
materials where the failure criterion is highly dependent on the mixed-mode ratio (as shown in 
Figure 3) the virtual crack closure technique (VCCT) [23, 24] has been widely used for 
computing energy release rates. Results based on continuum (2-D) and solid (3-D) finite element 
analyses provide the mode separation required when using the mixed-mode fracture criterion. 

The mode I, and mode II components of the strain energy release rate, Gi, Gu are 
computed as shown in Figure 4 for a 2-D four-noded element as an example of VCCT. The terms 
X’i , Z\ are the forces at the crack tip at nodal point i and u ’ e , w’ t (it ’ w ’,) are the displacements 
at the corresponding nodal points i and i behind the crack tip. For geometric nonlinear analysis 
where large deformations may occur, both forces and displacements obtained in the global 
coordinate system need to be transformed into a local coordinate system (x 1 , z 1 ) which originates 
at the crack tip as shown in Figure 4. The local crack tip system defines the tangential (x or 
mode II) and normal (z\ or mode I) coordinate directions at the crack tip in the deformed 
configuration. The extension to 3-D is straight forward and the total energy release rate GY is 
calculated from the individual mode components as GY =Gi +Gn +Gm, where Gm =0 for the 
two-dimensional case shown in Figure 4. 


2.2.2 A Global/Local Shell 3D Modeling Technique 

Built-up structures are traditionally modeled and analyzed using plate or shell finite 
elements as shown in Figure la to keep the modeling and computational effort affordable. 
Computed mixed mode strain energy release rate components, however, depend on many 
variables such as element order and shear deformation assumptions, kinematic constraints in the 
neighborhood of the delamination front, and continuity of material properties and section 
stiffness in the vicinity of the debond when delaminations or debonds are modeled with plate or 
shell finite elements [25-27]. These problems may be avoided by using three-dimensional 
models. Since many layers of brick elements through the thickness are often necessary to model 
the individual plies, however, the size of finite element models required for accurate analyses 
may become prohibitively large. 

For detailed modeling and analysis of the delaminations, the shell/3D modeling technique 
will reduce the modeling time since existing plate or shell models may be modified to shell/3D 
models. This is a considerable advantage compared to the creation of an entirely new three- 
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dimensional finite element model. The technique will also reduce computational time because 
only a relatively small section of interest needs to modeled with solid elements keeping the 
number of unknowns small. The technique combines the accuracy of the full three-dimensional 
solution with the computational efficiency of a plate or shell finite element model and has been 
demonstrated for various applications such as fracture toughness characterization specimens [28, 
29], on coupon level for the skin/stringer separation specimen [30, 31] and recently in a related 
study for skin/stringer separation [12] as shown in Figure 1. 


3. FINITE ELEMENT MODELING 

The stringer stiffened panel as shown in Figure 1 is bolted to a steel picture frame and 
subjected to pure shear loading. For the current study seven delamination lengths (a=8 1 .9 mm, 
88.9 mm, 94.9 mm, 101.6 mm, 203.2 mm, 279.4 mm and 355.6 mm) were selected to be 
reanalyzed with a smaller local 3D insert models to reduce the CPU time. Compared to the local 
insert models shown in Figure 1, only a small part of the intact and delaminated sections of the 
stringer were modeled with solid elements. 


3,1 Global Shell Model of Stringer Stiffened Panel 

The global model representing the steel loadframe, the panel made of graphite/epoxy 
prepreg tape and the stringers made of graphite/epoxy fabric was created using standard S4 shell 
elements available in the finite element software ABAQUS ' . The outer steel load frame was 
modeled with beam elements as shown in Figure 5a. The inner steel load frame which overlaps 
the panel edge was modeled with standard shell S4 elements. Shell elements were also used to 
model the panel and the stiffener components, such as foot, web and hat. Material data are given 
in Table 1. Details of the global finite element model are discussed in reference [12]. 

A detail of the global finite element model in the vicinity of the stringer termination is 
shown in Figure 5b. In preparation for the global/local modeling approach, shell elements 
representing the foot of the stiffener and the panel were removed from a small section of the 
original shell model around the center stringer tennination. The shell elements used to model the 
stiffener web and hat were kept in place. Shell edges (e.g. STRFOOT, STRBOT, PNLEND, 
PNLLEFT, PNLFRNT) in ABAQUS® were defined as shown, which were used to connect the 
shell model with the local 3D insert model. The connection was accomplished using the shell to 
solid coupling option in ABAQUS®, which allows the connection between non-conforming shell 
and solid models. In the detail it is also clearly visible that the panel skin and the stiffener foot 
are modeled as separate offset entities. The S4 shell elements are located at the panel skin and 
stiffener foot respective mid-planes. The shell elements are connected by beam elements 
designed to enforce plate theory constraints [12], In the sections containing the artificial defects 
the beam elements were replaced by gap elements. The gap elements allow the modeling of the 
skin/stringer separation but also prevent element interpenetration in case the surfaces get into 
contact during the analysis. 


3.2 Local 3D Insert Model for Solid Modeling of the Stringer Foot and Panel Skin 

The local 3D insert model was generated using C3D8I solid brick elements and consisted 
of an intact section and a delaminated section with a fine mesh around the delamination front as 
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shown in Figure 6a. Surfaces were defined in ABAQUS on the outer faces of the insert model to 
provide a connection with the global shell model using the shell to solid coupling option in 
ABAQUS” . The initial defect is located at the bondline between stringer foot and the panel. This 
defect was treated as a delamination and modeled as a discrete discontinuity using two 
unconnected nodes with identical coordinates one on each side of the delamination. A refined 
mesh was used along the stringer boundary in order to capture edge effects as shown in 
Figure 6b. Using the finite sliding option available in ABAQUS® contact was modeled between 
the delaminated surfaces to avoid interpenetration during analysis. Four elements over the 
thickness were used to model the foot of the stiffener made of carbon/epoxy fabric as shown in 
Figure 6b. The -45° skin ply made carbon/epoxy tape which is adjacent to the plane of 
delamination was modeled with one element. The remaining 10 plies of carbon/epoxy tape were 
modeled with three elements over the thickness a shown in Figure 6b. 


3.3 Combined Global/Local Shell/3D Model of Stringer Stiffned Panel 

For modeling the experiment, which was performed under constant displacement control, 
uniform displacements u,v were applied at one corner node to introduce shear as shown in 
Figure 7a. The in-plane displacements u,v were suppressed at the diagonally opposite corner and 
the out of plane displacements w were suppressed along all four edges across the entire width of 
the inner and outer steel load frame. The local 3D insert model containing a straight delamination 
front was inserted into the global shell model as shown in the detail of Figure 7b. The local 3D 
insert model consisted of an intact section and a delaminated section of length a with a fine mesh 
around the delamination front as discussed above. The global shell model was connected to the 
local 3D insert model using the shell-to-solid coupling option in ABAQUS® which allows the 
connection between non-conforming shell and solid models. For the entire analyses the non- 
linear solution option was used in ABAQUS 8 . 

A total of seven delamination lengths were modeled (a=81.9mm, 88.9 mm, 94.9 mm, 
101.6 mm, 203.2 mm, 279.4 mm and 355.6 mm). The initial length corresponds to the length of 
insert used to create an initial defect at the termination of the center stringer. Additional lengths 
were chosen to study the change in failure index distribution across the width ( b ) of the stringer 
with increasing delamination length (a). The lengths modeled correspond to the lengths 
investigated in a related study which was used for comparison [12]. A local 3D insert with a fine 
mesh near the free edges of the stringer foot is shown in Figures 8a. In the local model only 
38.1 mm of the stringer foot and panel skin were modeled with CD8I solid brick elements 
compared to 1 16.7 mm in the related study as shown in Figure lb. The number of element were 
reduced by 40% as listed in Table 2. Additionally a local 3D insert of the same length but with a 
uniform mesh across the entire width of the stringer foot was generated. The mesh was used for 
comparison and is shown in Figure 8b. The delaminated section of the stringer modeled with 
shell elements and the local 3D insert model used to analyze a delamination length <7=355.6 mm 
is shown in Figure 8c. In the local model only 55.6 mm of the stringer foot and panel skin were 
modeled with CD8I solid brick elements compared to 431.2 mm in the related study [12] as 
shown in Figure lc. The number of elements in the local insert could be reduced by nearly 60%. 
The size of the different solid insert models used is listed in Table 2 for comparison. 
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4. ANALYSIS RESULTS 

4.1 Deformed Finite Element Model for Solid Modeling of the Stringer Foot and Panel Skin 

The detail of the deformed local models are shown in Figures 9 and 10 for the last 
analysis increment after the entire external displacement u=v=6 . 35 mm had been applied. The 
detail of the deformed local models for a delamination length a=8 1.9mm are shown in Figure 9. 
In Figure 9a the opening gap elements are clearly visible. The gap elements were used to prevent 
contact in the delaminated section modeled with shell elements. For clarity the gap elements are 
not displayed for the remainder of the study. An opening of the delaminated section is observed 
across the width of the stringer over the entire delaminated section as shown in Figure 9b. The 
same deformation was also observed for the local insert model shown in Figures 9c where the 
stinger foot and adjacent panel skin were meshed uniformly with 30 elements across the width of 
the stringer foot. The same opening is also observed for modeled delamination lengths 
a= 88.9 mm and 94.9 mm as shown in Figures 10a and b. The detail of the deformed local model 
for a delamination length a= 1 0 1 .6 mm is shown in Figure 10c for the last analysis increment. 
The figure reveals that not the entire delaminated section opens under mode I. After initial 
opening, the section below the web termination closes and the delaminated surfaces contact. This 
closing is caused by a change in the local buckling pattern, due to stiffness changes caused by the 
longer delamination. It was observed that the local buckling pattern in the immediate 
surrounding of the delaminated stringer is very sensitive to the delamination length modeled, 
which made convergence difficult. For the models representing delamination lengths 
a=203.2 mm as show in Figure lOd and <2=279.4 mm as shown in Figure lOe an opening of the 
delaminated section is observed for the entire delamination length modeled. Partial closure of the 
delaminated surfaces which result in local contact is observed for the analysis of the <2=355.6 mm 
delamination as shown in Figure lOf. The observed deformations match the results from a related 
study where the entire delaminated section of the stringer had been modeled with solid elements 
[ 12 ]. 

4.2 Calculation of Mixed-Mode Strain Energy Release Rates and Failure Indices 

The virtual crack closure technique (VCCT) - discussed earlier - was used to calculate 
the mode contributions Gi, Gu and Gm the total energy release rate G j= Gi+Gn+Gm, as well as 
the mixed mode ratios Gs/Gt along the delamination front across the width b of the stringer for 
all delamination lengths modeled. Here Gs denotes the sum of the in-plane shearing components 
Gu+Giu . For two-dimensional analyses, where G///=0, this definition is equal to the commonly 
used definition of the mixed mode ratio, Gn / Gt . For three-dimensional analysis, which also 
yields results for the scissoring mode Gm, the modified definition of Gs is introduced since a 
mixed-mode failure criterion, which accounts for all three modes is currently not available. 

For each nodal point along the delamination front the critical energy release rate G c was 
calculated from the mixed mode failure criterion for IM7/8552 graphite/epoxy (Figure 3) 

( X 4.46 

for the computed mixed-mode ratio Gs/Gt at each point. Subsequently the failure index 
Gt/G c was determined with the assumption that delamination propagation occurs for 
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( 3 ) 



For all delamination lengths modeled, the computed failure indices were calculated for 
every fifth load increment plus the final increment and plotted versus the dimensionless 
coordinate s across the width of the stringer b 

= L; 0.0<5<1.0. (4) 

At the left edge of the stringer the nodal point coordinates are equal to y=yo which yields 
5=0.0 and the right edge nodal point coordinates are equal to y=yt, which results in 5=1.0 as 
depicted in Figure 7. 

For all delamination lengths modeled the calculated failure indices are shown in 
Figures 11 to 17. Since it was difficult to get corresponding results at the same load increment 
for different models, failure indices were calculated and plotted for the last analysis increment 
after the entire external displacement n=v=6. 35 mm had been applied. The calculated failure 
indices for delamination length a=81.9 mm are shown in Figure 11 for different models used. 
The distribution obtained from the inserts shown in Figure 9b and c is plotted together with the 
result obtained from the model of a related study [12] shown in Figure lb. The failure index 
peaked at the edges (5=0.0 and 5=1.0) with an additional peak around the center (5-0.5) 
underneath the stringer. The values in the center were higher for the model with 30 elements 
across the width. This can be explained since the mesh is finer in the center and is able to 
describe the failure index distribution in more detail. Overall, the results were in good agreement 
and showed that the computation of the failure index was insensitive to minor changes in the 
model. 

The failure indices computed for delamination length a=88.9, 94.6, 101.6, 203.2, 279.4 
and 355.6 mm are shown in Figures 12 to 17 for different models used. The distributions 
obtained from the inserts shown in Figure 10 are plotted together with the result obtained from 
the model of a related study [12], For delamination lengths a =88.9 mm and 94.6 mm the results 
were in good agreement as shown in Figures 12 and 13. 

For longer delaminations <3=101.6, 203.2, 279.4 and 355.6 mm the distributions across 
the width obtained from the short inserts follow the same trend as the results from the original 
study however differences in the center and near the edges become visible. The differences may 
be caused by the local shell to solid coupling. For the original models shown in Figures lb and c 
the local 3D insert models are connected to many shell nodes and closely follow the deformation 
of the global shell model. Away from the stinger termination the center stringer is modeled with 
a coarse shell mesh as shown in Figure 7a and therefore fewer shell nodes connect to the 
significantly shorter local insert as shown e.g. in Figure lOe. The deformation of the shorter 
insert may therefore not follow the deformed global shell model as well as the models which are 
connected to more shell nodes. 

For the models used for delamination length a=101.6 mm and 203.2 mm the mesh of the 
stringer web and foot and panel skin adjacent to the insert was manually refined as shown in 
Figures 10c and d. The modified mesh yielded failure indices across the width which were in 
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better agreement with the reference results [12] as shown in Figures 14 and 15. It is assumed that 
additional local refinement of the shell models near the flange edges would further improve the 
results. With respect to the distributions shown in Figures 16 and 17 local refinements of the 
stringer shell model are also expected to yield better agreements with the reference results [12] 
for delamination lengths a=279.4 and 355.6 mm. Additional analysis will be required to validate 
this assumption. 


5. SUMMARY AND CONCLUDING REMARKS 

The skin/stringer separation of a graphite/epoxy composite panel reinforced with three 
stringers and subjected to pure shear loading was studied using computational fracture analysis. 
The shear loading causes the panel to buckle and the resulting out-of-plane deformation initiates 
skin/stringer separation at the location of an embedded defect. The panel and surrounding load 
fixture were modeled with shell elements. 

A small section of the stringer foot and the panel in the vicinity of the delamination front 
were modeled with a local 3D solid model. A total of seven delamination lengths between 81.9 
and 355.6 mm were modeled. Across the width of the stringer foot the mixed-mode strain energy 
release rates were calculated using the virtual crack closure technique. A failure index was 
calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy 
material. Computed failure indices were compared to corresponding results from a related study 
where the entire delaminated section of the stiffener foot had been modeled with solid elements. 

Computed failure indices were in good agreement with results from models where the 
entire delaminated section of the stiffener foot had been modeled with solid elements. Changing 
the model, which had a coarse mesh in the center and a fine mesh near the edges to a unifonn 
mesh across the entire width of the stiffener foot, had a negligible influence on the results. The 
results showed that the failure indices were fairly insensitive to minor changes in the model. It 
was observed that the shell-to-solid connection influenced the computed failure indices. Local 
refinement of the shell model across the stringer foot and web improved the results. It is assumed 
that additional local refinement of the shell model near the flange edges would further improve 
the results. Additional analysis will be required to substantiate this assumption. 

Overall the study confirmed that the section modeled locally did not have to include the 
entire delaminated section. The use of a smaller local inserts reduced model size without 
compromising the computed failure indices. 
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Table 1. 


Material Properties 


IM7/8552 Unidirectional Graphite/Epoxy Prepreg 

£■11 = 150.0 GPa 

£22 = 10.7 GPa 

£33 = 10.7 GPa 

vi 2 = 0.33 

vi 3 = 0.33 

V 23 = 0.45 

Gi 2 = 4.8 GPa 

G 13 =4.8 GPa 

G23 = 3.4 GPa 

IM7/8552 Graphite/Epoxy Plain Weave Fabric 

£ll =73.1 GPa 

£22 = 73.1 GPa 

£33 = 10.6 GPa 

vi 2 = 0.04 

vi 3 = 0.35 

V 23 = 0.35 

Gi 2 = 4.8 GPa 

G 13 =4.1 GPa 

G 23 =4.1 GPa 

8552 Resin 

£=4.67 GPa 

v= 0.37 

(assumed isotropic) 

Steel 

£=206.8 GPa 

v= 0.33 

(assumed isotropic) 


The material properties are given with reference to the ply coordinate axes where index 1 1 denotes 
the ply principal axis that coincides with the direction of maximum in-plane Young’s modulus 
(fiber direction). Index 22 denotes the direction transverse to the fiber in the plane of the lamina 
and index 33 the direction perpendicular to the plane of the lamina. 
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Table 2. 


Size of Solid Finite Element Insert Model 


Delamination 
length a * in 
mm 

Local 3D insert, i 
length c * in mm 

-eference [12] 

Number of solid 
elements C3D8I 

Local 3 

length c * in mm 

D insert 

Number of solid 
elements C3D8I 

81.9 

Short local 3D 
insert 1 c=l 13.7 mm 

8320 

38.1 

4992 

88.9 

Short local 3D 
insert 1 c=l 13.7 mm 

7488 

39.2 

5200 

94.6 

Short local 3D 
insert 1 c=l 13.7 mm 

8112 

31.8 

4576 

101.6 

Short local 3D 
insert 1 c=l 13.7 mm 

8320 

23.8 

4160 

127.0 

Long local 3D 
insert 1 1 c=431.2 mm 

11440 

" 


203.2 

Long local 3D 
insert 1 1 c=431.2 mm 

12480 

63.5 

6656 

279.4 

Long local 3D 
insert 1 1 c=431.2 mm 

13728 

63.5 

6448 

355.6 

Long local 3D 
insert 1 1 c=43 1 .2 mm 

14768 

55.6 

6240 


*Delamination length a and length of local 3D insert c as defined in Figure 7b 
1 Short local 3D insert as shown in Figure lb 
1 Long local 3D insert as shown in Figure lc 
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a. Deformed three-stringer panel under shear loading 



b. Short local 3D insert near stringer termination with a=81.9 mm delamination length 



Figure 1: Stringer Stiffened Panel [12] 
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Interlaminar sliding shear 
Mode II 


Figure 2: F racture M odes. 



Mode III 
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PNLFRNT 


skin/stringer connection 
in intact section 
modeled with beam 
elements 


shell elements 
removed for local 
3D insert 


skin/stringer connection in 
delaminated section 
modeled with gap elements 


b. Detail of global shell model with edges for shell to solid coupling 
Figure 5: Global shell model of stringer stiffened panel 
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SOLTOP 



a. Local 3D insert model and surfaces for shell to solid coupling 


stringer foot made of 
IM7/8552 fabric [0/45/45/0] 
modeled with four solid 
elements 

skin made of IM7/8552 tape 
ply at interface modeled 
with one solid [-45] 


layered solid [-45/45/45] 


layered solid [0/0/0] 


layered solid [45/-45/-45/45] 



delamination front' 

delamination - 
delaminated section 


b. Detail of local 3D insert model around delamination front 


Figure 6: Loccd insert model 
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a. Finite element model with load and boundary conditions. 



b. Detail of center stringer with local 3D insert 
Figure 7: Global shellAocal 3D model of stiffened panel 
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a. Finite element model of local 3D insert for a=81.9 mm delamination length 




b. Local 3D insert for a=81.9 mm delamination length with uniform mesh 



c. Finite element model of local 3D insert for a=355.6 mm delamination length 
Figure 8: Finite element models for different delamination lengths 
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a. Deformed local 3D insert model with visible gap elements 


b. Deformed local 3D insert model with refined mesh at the edge 




c. Deformed local 3D insert with uniform mesh across the width 
Figure 9: Finite element model of local 3D insert for a=81.9 mm delamination length 
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c. Local 3D insert for a=101.6 mm d. Local 3D insert for a=203.2 mm 



e. Local 3D insert for a=279.4 mm f. Local 3D insert for a=355.6 mm 


Figure 10: Deformed center stringer with local 3D inserts for different delamination length a 
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O insert model (Fig. 1b) [12] 

• small insert (Fig. 9b) 

■ uniform 30 elements across width (Fig. 9c) 



0 0.2 0.4 0.6 0.8 1 


coordinate s, - 

Figure 11. Computed energy release rate across the width 
of the stringer obtained from different models for a=81.9 mm 
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0.4 0.6 

coordinate s, - 


Figure 14. Computed energy release rate across the width of the stringer for a= 101.6 mm 
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Figure 15. Computed energy release rate across the width of the stringer for a= 203.2 mm. 
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